Harmonic oscillations and their switching in 
elliptical optical waveguide arrays 



Ming Jie Zheng, 1 Yun San Chan, 1 and Kin Wah Yu 1 ' 2 

department of Physics, The Chinese University of Hong Kong, Shatin, New 
Territories, Hong Kong, China 

institute of Theoretical Physics, The Chinese University of Hong Kong, Shatin, New 
Territories, Hong Kong, China 

E-mail: m j zhengSphy . cuhk . edu . hk 

Abstract. We have studied harmonic oscillations in an elliptical optical waveguide 
array in which the coupling between neighboring waveguides is varied in accord with 
a Kac matrix so that the propagation constant eigenvalues can take equally spaced 
values. As a result, long-living Bloch oscillation (BO) and dipole oscillation (DO) are 
obtained when a linear gradient in the propagation constant is applied. Moreover, 
we achieve a switching from DO to BO or vice versa by ramping up the gradient 
profile. The various optical oscillations as well as their switching are investigated by 
field evolution analysis and confirmed by Hamiltonian optics. The equally spaced 
eigenvalues in the propagation constant allow viable applications in transmitting 
images, switching and routing of optical signals. 
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1. Introduction 

Optical oscillations in optical waveguide arrays (OWA) are harmonic when the 
propagation constants of the transverse modes are equally spaced [1, 2, 3]. These 
oscillator waveguide arrays are of particular interest not only due to their applications in 
transmitting images [4], focusing and steering [5], tunneling [6, 7, 8], and switching and 
routing optical signals [9, 10], but they are also good candidates for realizing the optical 
analogies of the dynamics in quantum systems [11]. Among these, Bloch oscillation (BO) 
and dipole oscillation (DO) are two fundamental and important types of oscillations 
[12, 13]. BO is known as the oscillatory motion of a particle in a periodic potential 
under a constant force [14], whose optical equivalent in an OWA is a linear gradient 
of the onsite propagation constants [15, 16]. DO usually occurs at the bottom of a 
parabolic band [17]. In the reciprocal space, DO distinguishes from BO that the BO 
momentum is accelerated monotonically by the force while the DO momentum oscillates 
due to the confinement of a band [18, 19]. The steering between BO and DO has been 
realized in a parabolic optical waveguide array [19]. However, both BO and DO decay 
fairly quickly due to the unequally spaced propagation constants in a parabolic band. 
In order to obtain long-living BO and DO, we need to construct OWAs so that the 
eigenvalues of the propagation constants are equally spaced. 

In the literature, several models have been proposed to achieve this goal, for 
example, the Wannier-Stark waveguide arrays, which require an infinitely varying 
transverse refractive-index distribution [15]. In real experiments [9, 10], we need 
to truncate the system, which unavoidably disrupt the equally spaced propagation 
constants. As a result, the oscillations can decay rapidly, especially when the light 
beam is incident at the boundaries of the system. Besides, there also exist dissipation 
[20], defects [21], and other kinds of damping effects in real systems. Recently, a finite 
waveguide array with graded coupling constants described by the elements of a Kac 
matrix was proposed [22, 23]. This is to achieve the required equally spaced propagation 
constants [3]. From now on, we refer these graded couplings k to circular couplings as k 
varies circularly as a function of waveguide index. In the finite waveguide arrays with a 
constant on-site potential, harmonic oscillation other than BO has been observed. The 
harmonic oscillation is similar to DO, as it also oscillates in the reciprocal space. 

In this work, we apply a linear gradient to the onsite propagation constants, which 
leads to the occurrence of BO. Strictly speaking BO is approximate, because the band 
structure is elliptical caused by the combination of linear graded onsite potential and 
the circular couplings. Due to the elliptical profile of the dispersion relation, we call 
the array as the elliptical optical waveguide array (EOWA). We propose to realize a 
switching between DO and BO or vice versa by ramping up the gradient of propagation 
constants. The DO-BO transition is investigated through the field-evolution analysis 
and confirmed by the Hamiltonian optics approach. The equally spaced propagation 
constants and switching between DO and BO have potential applications in optical 
steering devices. 
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2. Model and formalism 

The EOWA consists of N = 100 waveguides as shown schematically in figure 1. This 
array is divided into two zones (zone 1: < z < Z\ and zone 2: Z\ < z < z 2 ) along 
the longitudinal z-axis direction, where there are two elliptical potential profiles Hi and 
H 2 , respectively. As shown in figure 2, each elliptical potential profile is formed by 
the circular couplings between adjacent waveguides (k, solid line) and a linearly graded 
onsite propagation constant (f3 , dashed line). The switching between BO and DO can 
be achieved by ramping up the gradient of the propagation constant from zone 1 to zone 
2. The circular couplings are realized by carefully designing the structure of EOWA. 
The linear gradient in propagation constants are obtained by taking advantage of the 
electro-optical effect or thermo-optical effect [15, 9], and the the gradient strength is 
adjusted by the external voltage or temperature difference. The size of each waveguide 
is in the micrometer scale. However, the real physical parameters should be calibrated 
through experiments. As shown in figure 1, the intensity of the discrete input beam 
has a Gaussian distribution and its wave front is not a plane wave. The finite phase 
difference of excitations between adjacent waveguides is realized by putting a dielectric 
block just in front of the waveguide array. The light beam propagates along the axis 
of the waveguide array, that is, the z direction. The waveguide array is labeled by n 
(n — 1, 2, N) in the transverse direction. 

As mentioned in the literature [24], the waveguide arrays are analyzed by the 
coupled-mode theory [15]. The evolutionary equation of modal amplitude a n in the 
nth waveguide is written as 

{^lz + + Kn ^ ian ~ 1 ^ + K n,n+l a n+l(z) = , (1) 

where V n = an/N + Vo is the onsite propagation constant, a is the gradient of the linear 
potential, Vq is the onset propagation constant, which is set to be 4 in this study, and 
« n , n -i and K n , n+ i are the coupling constants between nearest-neighbor waveguides. In 
order to obtain harmonic BO and DO, the values of the coupling constants should be 
designed appropriately, as described in the following. 

Substituting the solution a™{z) = u™ exp(i(3 m z) into (1), we have 

0m< = Vn< + Kn,„-ltCi + «n,n+l<Vl , (2) 

where (3 m means the wavenumber of the supermode m and the transverse mode profile 
is given by a superposition of the mode amplitudes u™ of the individual waveguides. 
Equation (2) is rewritten in the matrix form 

0|u) = H|ii), (3) 

where the Hamiltonian matrix H is defined as H nn = V n = an/N + 4, H nn _x = « n;n _i, 
and H n>n+ i = K ntU+1 . In order to achieve the harmonic BO and DO, the eigenvalues 
of H should be equally spaced. It was shown that the iV + 1 by N + 1 Kac matrix 
K has equally spaced eigenvalues 2n — iV (n = 0, 1, ...,N) [22, 23]. The Kac matrix 
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is a tridiagonal matrix with elements K n+1 ^ n = n, K njJl+1 = N — n, and K nm = 
(otherwise), and is asymmetric (non-Hermitian). The Kac matrix is symmetrized [23] 
and the super and sub-diagonals of H are defined as [3] 

y/n(N - n) 

H n ,n+1 = H n+ i yn = K n ^ n+ i = K n+ i :n = — . (n — 1, 2, jV — 1)(4) 

The symmetric Hamiltonian matrix H is similar to the Kac matrix in a sense that 
both their eigenvalues are equally spaced. The column vector \u) and (3 denote the 
eigenvectors and eigenvalues of H, respectively. At a certain gradient a, the eigenvalues 
f3 have equal spacing A/3 = \fA + a 2 /N, and thus the density of states is a constant 
l/(NAp) = 1/V4 + CX 2 . 

Using the Hamiltonian matrix H, (1) is written as a z-dependent equation 

-iiu)=H|u). (5) 

dz 

It is analogous to the Schrodinger equation in quantum system, 

-i||0) = H|0). (6) 

Here h is taken to be unity. Thus the quantities (/3,z) in optical waveguide arrays 
corresponds to (u, t) in quantum system. We refer to the functional dependence of /3 on 
transverse wavenumber k as the dispersion relation in periodic optical waveguide arrays. 
For graded arrays, we can divide the infinite waveguide arrays into a large number of sub- 
waveguide arrays in the transverse direction, each of which can be regarded as infinite 
in size. Based on the treatment of graded system, we can obtain the band structure 
approximately as follows. The solution satisfies the relation u n+ \ = u n exp(ik), where 
k is the transverse wavenumber. Substituting this relation into (2), we obtain the 
approximate pseudo-dispersion relation 



f3(n,k,a)^a— + 4+ v y — J -cosk, (7) 

which is a function of waveguide index n, wavenumber k and gradient strength a. 

In the following calculations, we use the model parameters: the waveguide number 
iV = 100, the input waveguide number n = 40, the input wavenumber k = 0, the 
width of the input beam a = 5, the initial gradient strength a — 1.0 and the ramped-up 
gradient strength a = 4.0 in zone 1 and zone 2, respectively. 



3. Results 



3.1. Various oscillations and normal modes in EOWA 

To analyze the various oscillations in EOWA, we resort to an effective diagrammatic 
approach with the aid of a phase diagram [18]. From the pseudo-dispersion relation (7), 
two curves /3(n, 0) and /3(n, ir) serve as the upper and lower boundaries of the phase 
diagram for EOWA, as shown in figure 3(a). Separated by two critical lines /3 c i = (3(1, ir) 
and P C 2 = (3(N,0), there are three different regions: lower DO region, BO region, and 
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upper DO region. In these three regions, DO at the smaller n side, BO at the middle, 
and DO at the larger n side take place, respectively. As stated in the previous research 
[18], there exist correspondences between various optical oscillations and the localization 
of different gradon modes in OWA [18]. The mode patterns of all the eigenmodes are 
shown by the contour plots of square moduli of eigenmodes as a function of eigenvalue 
f3 and waveguide index n, as shown in Fig. 3(b). Separated by j3 c i and /3 c2 , there are also 
three different regions, which indicates three different gradon modes corresponding to 
the three kinds of optical oscillations. If we construct the light beam using components 
of a specific kind of gradon modes in a particular region, the light beam will undergo 
a certain type of oscillation, as shown by the three input Gaussian beams in different 
oscillation regions in figure 3(a). For a certain input beam, is it possible to undergo DO 
under certain conditions and BO under other conditions? The answer is yes, we will 
present it in detail in the following sections. 

3.2. DO-BO transition 

The switching between DO and BO is realized by ramping up the gradient of the 
propagation constant. Let us first sketch an example of the DO-BO transition as shown 
in figure 4(a). In the range < z < z±, the original onsite propagation constant has 
a linear gradient a = 1.0 and the Hamiltonian is Hi = k, 1.0) (solid lines). This 
linear gradient together with the circular couplings leads to the long-living DO for an 
input Gaussian beam (no = 40, k = 0, a = 5). In the original potential profile, the 
input light beam undergoes DO between points A and B. After two periods, the gradient 
of propagation constant is changed to be a = 4.0 when z\ < z < z 2 . As a consequence, 
the elliptical potential profile is changed to a new one (as the dashed lines showed), the 
light beam is lifted from point A to point C and undergoes BO between points C and 
D in the new potential profile. In a word, the DO-BO transition has been realized by 
varying the gradient strength of propagation constant from a — 1.0 to a — 4.0. As 
shown in figure 4(b), the DO-BO transition process A— > B — > C — > D is also marked on 
the phase-space orbit. The solid lines are for the original potential profile when a = 1.0 
and j3 = 5.38. Since these curves are closed and their wavenumbers are confined to a 
certain range, these features indicate the occurrence of DO. The dashed lines are for 
the new potential profile when a = 4.0 and /3 = 6.58. The value of wavenumber kd 
has no limitation along these dashed lines, which indicates the occurrence of BO. The 
transition takes place between points A and C, which coincide on the phase-space orbit. 

3.3. Field- evolution analysis 

The process of DO-BO transition is investigated through the field-evolution analysis. 
The analysis is performed with an input wave function at z — 0, 
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where k is the input transverse wave number. The incoming field at z (z < 0) is 
ip(z) = ip(0) exp(ip z), where (3q is the propagation constant of individual homogeneous 
channel. The intensity profile |?/>(0)| 2 has a discrete Gaussian distribution centered at 
the noth waveguide with spatial width a. This input beam is a discrete Gaussian beam, 
whose intensity distribution is schematically shown as the dashed line in figure 1. The 
exponential factor exp [— ik (n — n )] denotes the phase differences between input beams 
excited on the nth and the n th waveguides. Other than the Gaussian distribution, the 
input beam can take other functional forms. 

We expand the input wave function in terms of the supermodes \u m ), 

|^(0)> = J2 A m\Um), (9) 

m 

where A m = (u m \ip(0)) is the constituent component of the input Gaussian beam. The 
subsequent wave function at propagation distance z is given by 

\iP(z)) = J2 A me^ mZ \u m ). (10) 

m 

At a certain propagation distance z, the wave function in the reciprocal space can be 
obtained by taking the following Fourier transform 

WM)=JWn,*))]. (11) 
By using these wave functions, the mean value of n and k are obtained as 

(w> 



The field evolution of DO-BO transition is demonstrated by the contour plots of 
\ijj(n,z)\ 2 in the real space and \4>(k,z)\ 2 in the reciprocal space, as shown in Figs. 5(a) 
and 5(b), respectively. When < z < Zi, oscillatory motions take place in both real 
space and reciprocal space (in a confined range), which indicates the occurrence of DO 
in this range. When Zi < z < z 2 , there is still oscillatory motion in the real space (with a 
smaller amplitude), which in the reciprocal space, the wavenumber kd is accelerated and 
can cover the whole range [— 7r, it] in the reduced Brillouin zone. These features clearly 
show the occurrence of DO-BO transition. The period of these harmonic oscillations Z 
can be obtained from ip(z + Z) = ip(z) and (10), where the condition exp(iA/3Z) = 1 
should be satisfied. Thus the period is derived as Z = 2n/A/3 = 2irN/y/4 + a 2 . We 
need to emphasize that there is almost no leakage of field from the main orbit of field 
evolution. Both DO and BO can persist for a long propagation distance. This property 
is important in practical applications of propagating light in optical waveguide arrays. 



3.4- Hamiltonian optics 



The results of DO-BO transition obtained through the field-evolution analysis are now 
confirmed by the Hamiltonian optics approach. From the pseudo-dispersion relation 
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(7), the evolution of the mean position and the mean wavenumber of the beam can be 
solved by using the equation of motion 

dx d(3(x,k,a) dk d(3(x,k,a) 

dz dk dz dx 

where the auxiliary variable x = 2n/N — 1 (— 1 < rr < 1) is used to replace the discrete 
waveguide index n, and fl(x,k,a) = 4 + a(l + x)/2 + y/1 — x 2 cos k. The numerical 
Hamiltonian optics results of the mean position (x) and the mean wavenumber (kd) 
are shown by solid lines in Figs. 5(c) and 5(d), respectively The dashed lines and dots 
in Figs. 5(c) and 5(d) are the mean values of (x) and (kd) calculated from the field- 
evolution analysis results. The mean values obtained from two methods match very well 
with each other. 

4. Discussion and conclusion 

As suggested in Ref. [3], the experimental realization of the EOWA is similar to the 
experimental structures GaAs/AlGaAs [10]. To obtain the required n values, the spacing 
between the waveguides was carefully designed based on the evanescent codirectional 
couplings [25]. In the above investigation, we considered the nearest-neighbor couplings 
between waveguides. This consideration is valid when the field overlapping is significant 
between neighboring waveguides only. We can also extend this study to more general 
OWAs with higher-order couplings. As is shown in the reference [26], the second- 
order coupling in OWA leads to nontrivial BO. It is instructive to consider higher- 
order couplings and the effects of defects in EOWA. Isolated defects can always disrupt 
the equally spaced eigenvalues and lead to damping, but Longhi designed localized 
defects carefully by supersymmetric quantum mechanics and obtained undamped optical 
oscillations [21]. From the basic physics point of view, an intrinsic harmonic system 
can naturally help us to access the competing damping mechanisms, like defects and 
disorder or even nonlinearity. Thus in the future work, we will investigate the nature of 
the harmonic oscillation against these effects. 

In summary, we have studied the dipole oscillation and Bloch oscillation in 
the elliptical optical waveguide arrays, which possess a linear gradient in the onsite 
propagation constants and circular couplings between neighboring waveguides. The 
couplings between neighboring waveguides are designed according to the tridiagonal 
elements of a Kac matrix, which is to achieve equally spaced eigenvalues. As a 
consequence, both DO and BO are long-living, which is important to optical steering 
applications. We also proposed to realize the DO-BO transition or vice versa by 
ramping up the gradient of propagation constants. The spatial evolution of the DO-BO 
transition is demonstrated through the field-evolution analysis and confirmed by the 
Hamiltonian optics approach. The long-living DO and BO and their switching have 
viable applications in transmitting images, switching and routing of optical signals. 



Harmonic oscillations and their switching in elliptical optical waveguide arrays 8 
Acknowledgments 

This work was supported by RGC General Research Fund of the Hong Kong SAR 
Government. We thank Prof. K. Yakubo for careful reading of the manuscript and for 
many useful discussion and helpful suggestions. 



Harmonic oscillations and their switching in elliptical optical waveguide arrays 



9 



Schirmer S G, Fu H, and Solomon A I 2001 Complete controllability of quantum systems Phys. 
Rev. A 63, 063410 

Mirrahimi M and Rouchon P 2004 Controllability of Quantum Harmonic Oscillators IEEE Trans. 

Automat. Contr. 49, 745-747 
Gordon R 2004 Harmonic oscillation in a spatially finite array waveguide Opt. Lett. 29, 2752-4754 
Friesem A A, Levy U, and Silverberg Y 1983 Digital and analog optical broad-band transmission 

Proc. IEEE 71, 208 

Verslegers L, Catrysse P B, Yu Z, and Fan S 2009 Deep-subwavclength focusing and steering of 
light in an aperiodic metallic waveguide array Phys. Rev. Lett. 103, 033902 

Trompeter H, Pertsch T, Lederer F, Michaelis D, Streppel U, and Brauer A 2006 Visual 
Observation of Zcncr Tunneling Phys. Rev. Lett. 96, 023901 

Trompeter H, Krolikowski W, Neshev D N, Desyatnikov A S, Sukhorukov A A, Kivshar Y S, 
Pertsch T, Peschel U and Lederer F 2006 Bloch Oscillations and Zcncr Tunneling in Two- 
dimensional Photonic Lattices Phys. Rev. Lett. 96, 053903 

Dreisow F, Szameit A, Hcinrich M, Pertsch T, Nolte S, Tunnermann A, and Longhi S 2009 Bloch- 
Zener Oscillations in Binary Supcrlattices Phys. Rev. Lett. 102, 076802 

Pertsch T, Dannberg P, Elflein W, Brauer A, and Lederer F 1999 Optical Bloch oscillations in 
temperature tuned waveguide arrays Phys. Rev. Lett. 83, 4752-4755 

Morandotti R, Peschel U, Aitchison J S, Eisenberg H S, and Silberberg Y 1999 Experimental 
observation of linear and nonlinear optical Bloch oscillations Phys. Rev. Lett. 83, 4756-4759 

Longhi S 2009 Quantum-optical analogies using photonic structures Laser and Photon. Rev. 3, 
243-261 

Christodoulides D N, Lederer F, and Silberberg Y 2003 Discretizing light behaviour in linear and 

nonlinear waveguide lattices Nature 424, 817-823 
Lederer F, Stegemanb G I, Christodoulides D N, Assanto G, Scgev M, and Silberberg Y 2008 

Discrete solitons in optics Phys. Rep. 463, 1-126 
Hartmann T, Keck F, Korsch,H J and Mossmann S 2004 Dynamics of Bloch oscillations New J. 

Phys. 6, 2-25 

Peschel U, Pertsch T, and Lederer F 1998 Optical Bloch oscillations in waveguide arrays Opt. 
Lett. 23, 1701-1703 

Sapicnza R, Costantino P, Wiersma D, Ghulinyan M, Oton C J, Pavesi L 2003 Optical analogue 

of electronic Bloch oscillations Phys. Rev. Lett. 91, 263902 
Menotti C, Smerzi A, and Trombettoni A 2003 Supcrfluid dynamics of a Bose-Einstcin condensate 

in a periodic potential New J. Phys. 5, 112.1-112.20 
Zheng M J, Xiao J J, and Yu K W 2010 Controllable optical Bloch oscillation in planar graded 

optical waveguide arrays Phys. Rev. A 81, 033829 
Zheng M J, Chan Y S, and Yu K W 2010 Steering between Bloch oscillation and dipolc oscillation 

in parabolic optical waveguide arrays J. Opt. Soc. Am. B 27, 1299-1304 
Efremidis N K and Christodoulides D N 2004 Bloch oscillations in optical dissipative lattices Opt. 

Lett. 29, 2485-2487 

Longhi S Bloch oscillations in tight-binding lattices with defects 2010 Phys. Rev. B 81, 195118 
Taussky O and Todd J 1991 Another look at a matrix of Mark Kac Linear Algebra Appl. 150, 
341-360 

Edclman A and Kostlan E 1994 The road from Kac's matrix to Kac's random polynomials 
Proceedings of the Fifth SIAM on Applied Linear Algebra J G. Lewis, Ed, Philadelphia, 503 

Trompeter H 2006 Discrete optics in inhomogeneous waveguide arrays, PhD dissertation, 
Fricdrich-Schiller Universitat Jena 

Yariv A 1997 Optical Electronics in Modern Communications, Oxford University, New York, 5th 
edition 

Wang G, Huang J P, and Yu K W 2010 Nontrivial Bloch oscillations in waveguide arrays with 
second-order coupling Opt. Lett. 35, 1908-1910 



Harmonic oscillations and their switching in elliptical optical waveguide arrays 



10 



List of Figure Captions 

Fig. 1 (Color online) Schematic diagram for the elliptical optical waveguide arrays and 
the discrete input beam. The light propagates along the axis of waveguide, that is, the 
z direction. The waveguide array is labeled by n (n = 1,2, ..,N). The amplitude of 
the input beam has a Gaussian distribution. The phase differences between different 
waveguides are finite, which can be obtained by putting a dielectric block just in front 
of the array. EOWA.eps 

Fig. 2 (Color online) In the elliptical optical waveguide arrays, the elliptical potential 
profile is formed by the circular couplings (solid line) and a linear onsite propagation 
constant (dashed line). Coup.eps 

Fig. 3 (Color online) (a) Phase diagram for the elliptical optical waveguide array with 
N = 100 waveguides. The pseudo-dispersion relation lines /3(n, 0) and /3(n,ir) serve as 
the upper and lower boundaries of the band, whose shape is elliptical. Separated by 
the two critical lines f3 c i = /?(1, 7r) and /3 c2 = /3(n, 0), three regions are formed, which 
indicate three different kinds of oscillations: lower dipole oscillation, Bloch oscillation 
and upper dipole oscillation, (b) The contour plots of mode patterns of the eigenvectors 
in the real space. The shape of the contour plots is similar to the elliptical band. The 
corresponding mode patterns for the three different oscillations are also separated by 
the two critical lines (5 C \ and f3 C 2- PDMP.eps 

Fig. 4 (Color online) (a) A possible DO-BO transition realized by ramping up the value 
of a from 1.0 to 4.0. In the original potential profile, the light beam undergoes DO 
between points A and B. After several periods, the potential profile is lifted, the light 
beam is transited from point A to point C and undergoes BO between points C and 
D. (b) The phase-space orbits for the elliptical optical waveguide arrays when a — 1.0 
(solid lines) and a = 4.0 (dashed lines). The DO-BO transition process A— > B — > C — > 
D points A and B are also shown. DO-BO. eps 

Fig. 5 (Color online) Contour plots of field-evolution analysis results for (a) \ip(n, z)\ 2 as 
a function of the waveguide index n and the propagation distance z and (b)\(f>(k, z)\ 2 as 
a function of the transverse wavenumber k and the propagation distance z. Comparison 
of Hamiltonian optics results with field-evolution analysis results for (c) (x) and (d) 
(kd) in the DO-BO transition. Note that the propagation distance is rescaled by N/2. 
xkComp.eps 



Harmonic oscillations and their switching in elliptical optical waveguide arrays 



11 




Figure 1. (Color online) Schematic diagram for the elliptical optical waveguide arrays 
and the discrete input beam. The light propagates along the axis of waveguide, that is, 
the z direction. The waveguide array is labeled by n (n = 1, 2, ... , N). The amplitude 
of the input beam has a Gaussian distribution. The phase differences between different 
waveguides are finite, which can be obtained by putting a dielectric block just in front 
of the array. EOWA.eps 
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Figure 2. (Color online) In the elliptical optical waveguide arrays, the elliptical 
potential profile is formed by the circular couplings (solid line) and a linear onsite 
propagation constant (dashed line). Coup.eps 
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Figure 3. (Color online) (a) Phase diagram for the elliptical optical waveguide 
array with N — 100 waveguides. The pseudo-dispersion relation lines /?(«, 0) and 
/3(n, it) serve as the upper and lower boundaries of the band, whose shape is elliptical. 
Separated by the two critical lines f3 c i = /3(1,tt) and (3 C 2 = /3(n, 0), three regions are 
formed, which indicate three different kinds of oscillations: lower dipole oscillation, 
Bloch oscillation and upper dipole oscillation, (b) The contour plots of mode patterns 
of the eigenvectors in the real space. The shape of the contour plots is similar to the 
elliptical band. The corresponding mode patterns for the three different oscillations 
are also separated by the two critical lines /3 c i and /3 C 2- PDMP.eps 
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Figure 4. (Color online) (a) A possible DO-BO transition realized by ramping up the 
value of a from 1.0 to 4.0. In the original potential profile, the light beam undergoes 
DO between points A and B. After several periods, the potential profile is lifted, the 
light beam is transited from point A to point C and undergoes BO between points 
C and D. (b) The phase-space orbits for the elliptical optical waveguide arrays when 
a = 1.0 (solid lines) and a = 4.0 (dashed lines). The DO-BO transition process A— >• 
B — > C — >• D points A and B are also shown. DO-BO. eps 
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Figure 5. (Color online) Contour plots of field-evolution analysis results for (a) 
\ip(n, z)\ 2 as a function of the waveguide index n and the propagation distance z and 

(b) \(f>(k, z)\ 2 as a function of the transverse wavenumber k and the propagation distance 
z. Comparison of Hamiltonian optics results with field-evolution analysis results for 

(c) (x) and (d) (kd) in the DO-BO transition. Note that the propagation distance is 
rescaled by N/2. xkComp.eps 



